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SUMMARY 

With  the  objective  of  developing  analytical  traction  models  for  the  simula¬ 
tion  of  rheological  behavior  of  lubricating  oils  used  In  practical 
applications,  such  as  rolling  bearings,  gears,  cams  and  similar  mechanical 
components,  the  feasibility  of  a  traction  model,  based  on  viscosity- 
pressure-temperature  relations  in  elastohydrodynamic  contacts,  is 
demonstrated  for  two  military  lubricants  with  the  MIL-L-23699  and  MIL-L-7808 
specifications.  It  is  shown  that  the  model  predictions  agree  fairly  well 
with  the  experimental  data,  and  the  practical  utility  of  the  model  is 
demonstrated  by  incorporating  the  traction  model  in  the  advanced  bearing 
dynamics  computer  model,  ADORE,  and  by  evaluating  the  roller  slip  and  skid 
behavior  in  a  cylindrical  roller  bearing  as  a  function  of  the  lubricant 
traction  parameters. 

Two  types  of  rheological  models  are  considered:  type  I  where  the  viscosity 
varies  exponentially  with  pressure  and  temperature,  and  type  II  where  the 
variation  in  viscosity  is  exponential  with  pressure  and  inverse  temperature. 
With  the  prescribed  general  form  of  the  viscosity  behavior  the  energy  equa¬ 
tion  is  solved  through  the  lubricant  film,  with  appropriate  boundary 
conditions,  and  the  shear  stress  is  integrated  to  compute  the  total  traction 
in  the  contact.  The  computed  traction  is  compared  to  the  experimental  data 
and  a  nonlinear  regression  analysis  is  carried  out  to  determine  the  model 
coefficients  by  minimizing  the  squared  deviation  between  model  predictions 
and  experimental  observations. 

Experimental  data  obtained  at  the  AFWAL  Aero  Propulsion  Laboratory  for  two 
types  of  contacts  is  considered  for  determining  the  model  coefficients:  an 
elliptical  contact  with  a  set  of  cylindrical  disk  specimens  with  a  rolling 
radius  of  38.10  mm  and  a  crown  radius  of  914.^0  mm,  and  a  circular  contact 
with  a  set  of  disk  with  spherical  surfaces  with  a  radius  of  38.10  mm.  The 
lubricant  temperature  at  inlet  varied  from  2 94°  to  366  °K,  the  rolling 
velocity  ranged  between  21  and  M6  M/Sec,  and  the  Hertz  contact  pressures 
varied  between  500  and  2000  MPa.  The  elliptical  contacts  represent  the  low 
end  while  the  circular  contacts  correspond  to  the  high  end  of  Hertz  contact 
pressures. 


Both  the  type  I  and  type  II  viscosity  relations  fit  the  experimental 
data  obtained  with  the  elliptical  contact  configuration  for  both  the  MIL-L- 
23699  and  MIL-L-7808-type  lubricants  fairly  well.  For  the  type  I  model  the 
regression  coefficients  indicate  that  the  entire  variation  may  be  repre¬ 
sented  by  a  variation  in  a  reference  viscosity  parameter;  thus  a  greatly 
simplified  relation  for  traction  may  be  obtained.  Also,  since  the  the  math¬ 
ematical  simplicity  of  the  type  I  relation  allows  for  an  effective 
incorporation  of  the  traction  model  in  sophisticated  component  models,  such 
as  the  advanced  bearing  dynamics  computer  code  ADORE,  it  offers  significant 
potential  in  the  design  and  performance  optimization  of  advanced  mechanical 
components. 


In  the  case  of  circular  contacts  with  relatively  high  pressures,  there  is 

noticeable  scatter  between  the  model  predictions  and  experimental  observa-  - - 

tions,  particularly  in  terms  of  the  very  high  traction  slope  in  the  very  low  oc*es 


or 


slip  region.  Further  model  refinements  to  allow  for  the  viscoelastic  be¬ 
havior  of  the  lubricant  are  required  for  a  more  realistic  modeling  of  the 
lubricants  under  such  high  pressures. 

The  practical  significance  of  the  traction  model  is  demonstrated  by  carrying 
out  a  simulation  of  roller  slip  or  skid  in  a  cylindrical  roller  bearing  by 
incorporating  the  traction  model  in  the  bearing  dynamics  computer  program, 
ADORE.  The  very  low  traction  coefficients  with  the  MIL-L-7808  lubricant 
result  in  substantial  slip  of  the  roller  as  it  travels  through  the  contact 
zone  in  a  radially  loaded  bearing.  The  viscous  heating  in  the  lubricant 
film,  resulting  from  excessive  slip,  results  in  a  significant  reduction  in 
lubricant  film  thickness,  and  a  metal  contact  is  eminent  under  the  maximum 
load  point.  Such  a  close  coupling  between  the  bearing  performance  and  trac¬ 
tion  behavior  of  the  lubricant  is  further  demonstrated  by  simulating  a 
hypothetical  lubricant  with  somewhat  Increased  traction  and  by  showing  that 
the  roller  slip  behavior  may  be  significantly  improved  by  a  careful  altera¬ 
tion  of  the  lubricant  traction  behavior.  Thus  the  feasibility  of  developing 
the  optimum  lubricant  behavior  for  advanced  applications,  through  traction 
modeling  and  dynamic  simulation  of  component  behavior,  is  demonstrated. 
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INTRODUCTION 


The  significance  of  lubricant  behavior  in  the  performance  of  Intricate 
mechanical  components,  such  as,  rolling  bearings,  gears  and  transmis¬ 
sions,  has  been  well  established.  It  is  well-known  that  proper  lubrication 
prevents  excessive  metallic  contacts  and  reduces  friction  and  wear  of 
interacting  elements,  and  serves  as  an  effective  cooling  media.  The 
traction  behavior  of  the  lubricant  has  proven  to  be  a  key  factor  in 
determining  the  stability  and  the  overall  dynamic  performance  of  a 
mechanical  system.  For  example,  in  rolling  bearings  the  tractive 
forces  at  the  rolling  element/race  interface  greatly  influence  the  orbi¬ 
tal  acceleration  of  rolling  elements  which,  in  turn,  determines  the 
extent  of  rolling  element/cage  collisions  and  associated  cage 
instabilities.  Thus,  cage  instabilities,  ball  and  roller  skid,  roller 
skew  and  excessive  wear  associated  with  the  highly  dynamic  interaction 
between  the  bearing  elements  are  often  attributed  to  the  traction  behavior 
of  the  lubricant. 

The  need  for  a  clear  understanding  of  the  rheological  behavior  of  the 
lubricant,  and  the  development  of  realistic  models  which  can  accurately 
predict  the  lubricart  behavior  for  a  prescribed  application,  is  further 
demonstrated  by  apparently  conflicting  requirements.  For  example,  in 
order  to  minimize  direct  metal  contact  and  wear,  the  lubricant  film  thick¬ 
ness  between  the  interacting  elements  should  be  maximized,  which  requires 
the  lubricant  to  have  a  relatively  high  viscosity;  but  the  requirements  of 
low  traction,  heat  generation  and  torque  lead  to  a  low-viscosity  lubricant. 
In  order  to  simultaneously  satisfy  such  performance  requirements,  an 
intricate  coupling  between  the  lubricant  behavior  and  the  performance  of 
a  machine  element  under  consideration  is  necessary.  In  other  words,  the 
lubricant  model  must  be  incorporated  in  the  component  model  which  simu¬ 
lates  the  overall  performance  of  the  mechanical  component,  such  as  a 
rolling  bearing.  From  a  practical  standpoint,  this  leads  to  a  dual  sig¬ 
nificance  of  lubricant  traction  modeling:  first,  the  model  can  provide 
realistic  designs  of  machine  components;  and  second,  the  parametric  perfor¬ 
mance  evaluations,  which  can  be  obtained  by  the  combined  lubricant  and 
component  models,  can  help  identify  crucial  lubricant  properties  re¬ 
quired  for  a  given  performance  specification.  Thus  for  the  military 
applications  of  rolling  bearings,  where  the  stability  of  the  bearing 
elements  and  an  acceptable  performance  of  the  bearing  has  to  be  ensured 
under  the  most  adverse  operating  conditions,  a  realistic  modeling  of  the 
traction  behavior  of  the  lubricant  is  extremely  important. 
Furthermore,  the  required  understanding  of  lubricant  behavior  and  the 
development  of  predictive  models  Is  vital  to  future  DOD  needs  that  are 
beyond  the  current  technological  limits. 

Over  the  past  three  decades  a  significant  fraction  of  the  mechani¬ 
cal  engineering  literature  has  been  devoted  to  elastohydrodynamic 
lubrication  and  lubricant  rheology.  Investigations  have  ranged  from  the 
most  fundamental  problem  of  lubricant  behavior  at  a  molecular  level  to 
very  applied  problems  where  the  global  hydrodynamics  intricately  couples 
with  the  elastic  distortion  of  the  load-bearing  contacts  in  a  rolling 
bearing.  The  rolling-disk  type  of  test  rigs,  such  as  those  at  the  AFWAL 
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Aero  Propulsion  (AFWAL/POSL)  and  Materials  (AFWAL/MLBT)  laboratories,  have 
proven  to  effectively  simulate  a  concentrated  contact  and  to  provide  a 
means  for  evaluating  the  tractive  behavior  of  a  lubricant.  With  the  use 
of  such  test  rigs,  traction  data  for  a  number  of  military  lubricants  have 
been  obtained  [1-4]*.  The  commonly  used  military  oil,  with  MIL-L-7808 
specification,  has  been  tested  by  a  number  of  investigators.  However,  due 
to  both  the  complexities  in  lubricant  behavior  and  the  experimental  dif¬ 
ficulties  in  traction  measurement,  the  range  of  operating  conditions 
over  which  the  data  have  been  obtained  is  limited  and  the  repeatability  of 
the  data  has  often  been  questionable.  In  view  of  these  difficulties  the 
development  of  a  realistic  model  has  been  extremely  difficult,  and  the 
available  models  [1,3,5-10,25,26]  need  further  experimental  validation 
before  they  can  be  adopted  for  military  applications.  Also,  the  analytical 
formulations  require  substantial  refinement  before  the  models  can  be  effec¬ 
tively  used  in  bearing  dynamics  computer  models.  With  significant 
improvements  in  the  experimental  techniques,  the  recent  data  obtained 
at  the  AFWAL  Aero  Propulsion  Laboratory  for  the  MIL-L-7808-type  lubricant 
probably  provides  the  most  extended  and  reliable  data  base  for  the 
development  of  traction  model  for  such  military  lubricants. 

Correlations  of  the  traction  data  to  the  fundamental  properties  of  the 
lubricant,  requires  that  the  rheological  behavior  be  measured  at  very 
high  pressures,  typical  of  a  concentrated  contact  in  a  rolling  bearing. 
This,  also,  is  an  extremely  difficult  task,  and  the  available  lubricant 
rheology  data  at  high  pressures  is  very  scarce.  Again,  the  facilities 
at  the  AFWAL  Aero  Propulsion  Laboratory,  such  as  the  high-pressure 
viscometer,  indeed,  offer  the  most  unique  and  advanced  experimental 
apparatus  for  investigating  the  fundamental  behavior  of  lubricants  in 
concentrated  contacts. 

Along  with  the  advances  in  lubrication  theories,  the  computer  modeling 
of  rolling  bearing  dynamics  has  significantly  advanced  over  the  past  few 
years.  The  computer  models  BASDAP  [11],  BDYN  [12],  DREB  [13,14], 
RAPIDREB  [15-17],  TRIB01  [18]  and  more  recently  ADORE  [19]  have  become  well 
known,  and  their  capability  in  modeling  truly  dynamic  behavior  of  the 
bearing  under  the  most  sophisticated  and  adverse  operating  conditions  has 
been  greatly  demonstrated  [19-22].  For  an  effective  simulation  of  bearing 
behavior,  the  lubricant  traction  model  has  been  identified  as  a  key 
input  to  these  bearing  dynamics  computer  models  [14,17,19,21].  This 
further  emphasizes  the  need  for  the  development  of  realistic  traction  models 
for  various  lubricants. 

With  particular  emphasis  on  advanced  DOD  applications  of  rolling  bea¬ 
rings,  the  practical  significance  of  a  traction  model  must  be 
evaluated  in  terms  of  its  utility  in  the  bearing  dynamics  computer 
models.  Due  to  the  complexities  associated  with  the  dynamic  interaction 
between  bearing  elements,  the  computer  models  which  simulate  the  real-time 
dynamic  performance  of  a  bearing,  such  as  the  ones  listed  above,  are 
quite  complex,  and  they  very  often  require  a  substantial  computational 
effort.  Since  the  lubricant  model  operates  in  the  inner-most  loop  of 
these  computer  codes,  a  slight  computational  complexity  in  the  lubricant 
model  greatly  amplifies  the  overall  computational  effort,  and  it  imposes  a 
severe  restriction  on  the  effective  and  economical  use  of  both  the  bear¬ 
ing  dynamics  and  the  lubricant  traction  models.  It  is,  therefore, 
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necessary  that  the  lubricant  behavior  be  modeled  such  that,  along  with  a 
very  realistic  prediction  of  traction,  the  model  offers  a  high  computa¬ 
tional  efficiency.  It  is  this  practical  requirement  which  leads  to  a 
rather  limited  utility  of  the  available  traction  models,  and  the  need  for 
further  development  becomes  quite  clear. 

The  general  approach  to  modeling  the  lubricant  behavior  has  been  a 
two-step  process.  The  traction  behavior  is  first  experimentally 
measured  over  the  expected  range  of  operating  conditions  and  then  "curve 
fitted"  to  a  mathematical  model.  Some  of  the  models  may  essentially 
provide  an  interpolation  of  the  actual  data  while  the  others  may  lead  to 
certain  constitutive  equations  which  have  a  physical  significance.  The 
ultimate  requirement  is  that  the  model  be  computationally  "simple"  and  yet 
provide  a  realistic  prediction  of  the  traction  behavior.  With  such 
guidelines  the  objective  of  the  present  program  is  to  develop  effective 
traction  models  for  military  oils,  with  the  available  traction  data  as  a 
starting  point. 
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ANALYTICAL  APPROACH 


The  mechanics  of  an  elastohydrodynamic  contact  is  schematically 
described  in  figure  1.  Since  the  size  of  the  contact  is  generally  quite 
small  compared  to  the  radius  of  the  interacting  elements,  the  assumption  of 
a  relatively  flat  contact,  as  shown  in  figure  1,  is  quite  realistic.  There 
are  essentially  three  zones  in  the  contact  region:  inlet  zone,  high- 
pressure  contact  zone  and  the  exit  zone.  The  lubricant  flows  from  the 
inlet  to  exit  through  the  high-pressure  contact.  The  properties  of  the 
lubricant  in  the  inlet  zone,  which  is  essentially  at  ambient  pres¬ 
sure,  determine  the  nominal  thickness  of  the  lubricant  film  in  the 
contact,  while  the  tractive  behavior  is  determined  by  the  lubricant 
properties  under  the  high  pressure  exerted  in  the  contact.  In  the  exit 
zone,  phenomena,  such  as  cavitation,  are  significant,  but  for  the  computa¬ 
tion  of  lubricant  film  thickness  and  traction  it  is  reasonable  to 
restrict  the  modeling  to  the  inlet  and  high-pressure  contact  zones. 

A  realistic  simulation  of  traction,  therefore,  demands  that  the  lubricant 
properties  be  known  both  under  ambient  and  high-pres3ure  conditions. 
Although  the  viscosity-temperature  relations  under  ambient  or  very  low  pres¬ 
sures  can  be  easily  measured,  high-pressure  viscometry  still  requires 
considerable  development;  perhaps,  the  high-pressure  viscometer  at  the 
AFWAL  Aero  Propulsion  Laboratory,  will  provide  such  data  in  the  near 
future.  For  the  present  it  may  be  reasonable  to  assume  a  certain  type  of 
viscosity-pressure-temperature  relation  in  the  high-pressure  contact  and 
calculate  the  traction  as  a  function  of  the  coefficients  in  the  constitu¬ 
tive  relation.  The  computed  traction  can  be  compared  to  the  experimental 
data,  and  the  coefficients  which  provide  a  reasonable  fit  can  be  determined 
by  carrying  out  a  regression  analysis.  Such  an  approach  is  very  complemen¬ 
tary  to  the  high-pressure  viscometry  research.  The  computed  coefficients, 
when  compared  with  those  derived  for  the  high-pressure  viscometry  data,  as 
and  when  it  becomes  available,  shall  provide  significant  insight  into  the 
fundamental  behavior  of  the  lubricant  in  a  concentrated  rolling-sliding 
contact. 

The  general  approach  to  traction  model  development  in  the  present  investiga¬ 
tion,  therefore,  consists  of  the  following  steps: 


1.  Computation  of  lubricant  film  thickness. 

2.  Computation  of  traction. 

3.  Regression  analysis. 


2.1  Computation  of  Lubricant  Film  Thickness 

The  lubricant  film  thickness,  h,  in  an  elastohydrodynamic  contact  is 
generally  defined  as 
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Figure  1.  Schematic  Representation  of  an  Elastohydro- 
dynamic  Contact  Zone 
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with  a  -  Semi-major  axis  of  the  contact  ellipse 

b  -  Semi-minor  axis  of  the  contact  ellipse 

E  -  Elastic  modulus  of  the  interacting  body 

Q  -  Applied  load 

R  -  Radius  in  the  rolling  direction 
a  -  Pressure-Viscosity  coefficient 
v  -  Poisson’s  ratio 
M,  -  Absolute  viscosity 

and  the  subscripts  ,  and  2  denote  the  two  elastic  bodies  in  contact. 

The  thermal  reduction  factor,  $T,  may  be  determined  from  the  numerical 

results  tabulated  by  Smith  et  al  CIO],  based  on  the  work  by  Cheng  and 
Sternlicht  [24], 
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2.2  Computation  of  Traction 

The  computation  of  traction  is  based  on  a  model  originally  presented  by 
Kannel  and  Walowit  [53  and  later  used  by  Gupta  et  al.  [33.  The  problem  is 
basically  formulated  in  terms  of  the  following  equations: 


Energy  Equation: 

„3aT 

K3p“ 

• 

■  -  “T  3 

(3) 

Geometric  Compatibility: 

3u 

3z 

-  s(t,p,T) 

(«) 

Constitutive  Equation: 

s(t, 

P,T)  “  y(p,T) 

(5* 

where 

K  -  thermal  conductivity 
T  -  absolute  temperature 
t  -  shear  stress 

s  -  shear  strain  rate 
u  -  lubricant  velocity 
p  -  pressure 
y  -  lubricant  viscosity 

z  -  coordinate  across  the  film,  as  shown  in  figure  1. 

With  the  nominal  film  thickness  computed  by  equations  (1)  and  (2),  the  ap¬ 
plicable  boundary  conditions,  for  the  solution  of  equations  (3)  to  (5)  are 


z  *  0:  T  -  Tt  and  u  -  U, 

z  -  h:  T  -  T2  and  u  -  U2 


(6) 


Let  the  temperatures  of  the  two  contacting  bodies  be  equal,  i.e.,  T,-T2-Te, 
and  with  u0  as  the  viscosity  at  reference  temperature  T0,  all  the  variables 
may  be  nondimensionalized  as  follows: 


T*  -  T/T0,  y*  -  y/y0 i  z*  -  z/h  and  U»  -  u/U  (7) 


Equations  (3~5)  may  now  be  combined  to  obtain  the  following  governing 
equations: 


3*T  tUh  3U 

JtJ  KT0  3z 

and 

3U  _  hT  1 
3z  ’  y0U  y 


(8) 

(9) 


where  the  asterisk  has  been  dropped  for  brevity.  Also  note  that  y  is  a 
function  of  pressure  and  temperature  as  defined  in  equation  (5). 


Integrating  equation  (8)  with  the  first  of  boundary  condition  in  (6)  gives 
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(10) 


3T  _  3T 
3z  3z 


z-0 


i]7h 

KT 


-  U.) 


and  the  division  of  equation  (10)  by  (9)  results  in  the  following  governing 
equation 


1  II 

y  3U 


Hal?  3T 
hi  3z 


(U  -  U,) 


(11) 


which  may  be  Integrated,  through  the  lubricant  film,  with  the  boundary  con¬ 
ditions  in  equation  (6)  to  obtain 


3T 

3z 


z-0 


hiU(U»  -  Ut) 
2KT  a 


(12) 


By  substituting  equation  (12)  in  (11)  and  by  denoting  the  viscosity  as  a 
function  of  pressure  and  temperature 

y  -  f (p,T)  (14) 

the  governing  equation  through  the  lubricant  film  may  be  written  as 

i  -  -ftr <>  -  u>  r<p-T>  (,5> 

Also,  the  shear  stress  through  the  film,  is  obtained  by  combining  equation 
(4)  and  (5) 


fUa 

f(p,T)  dU  (16) 

Equations  (15)  and  (16)  are  solved  simultaneously  to  obtain  the  shear  stress 
through  the  film.  Equation  (15)  is  first  Integrated  to  obtain  a 
temperature-velocity  relation,  which  when  combined  with  (14)  gives  the 
viscosity-velocity  relation  required  in  (16). 

For  any  arbitrary  viscosity-pressure-temperature  relation,  equation  (14), 
equations  (15)  and  (16)  may  be,  in  general,  solved  numerically  to  obtain  the 
shear  stress  at  any  point  in  the  contact  zone.  The  shear  stress  can  then  be 
integrated  over  the  contact  zone  to  obtain  the  total  tractive  force.  A 
traction  coefficient  may  then  be  defined,  either  in  terms  of  the  local 
stress  ratio  or  the  ratio  of  the  total  traction  to  applied  normal  force 

tc,  ,  -  x/p 

local  * 

(17) 

tc  ,  .  .  -  F/Q 
global 

where  F  is  the  integrated  traction  force  and  Q  is  the  applied  normal  contact 
load. 
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Since  the  slip  velocity  varies  greatly  in  a  typical  rolling-element-to-race 
contact,  the  local  traction  coefficient  is  most  suitable  when  applying  the 
traction  model  to  a  rolling  bearing.  However,  for  the  purpose  of  laboratory 

investigations  to  determine  the  model  coefficients  with  single  contacts,  the 

global  traction  coefficient  is  significant. 

It  is  clear  from  the  above  formulation  that  equation  (1M),  which  defines  the 
pressure-viscosity  relation  of  the  lubricant,  is  the  primary  input  to  the 
traction  model.  The  complexities  associated  with  traction  computation  are 
directly  related  to  the  general  form  of  this  constitutive  equation.  The  two 
most  common  forms  of  the  pressure-temperature-viscosity  relation  are: 

Type  I  relation:  -4r  -  f(p,T)  -  exp  {a*p  +  B*T0 ( 1 -T) }  (18) 

Uo 

Type  II  relation:  — *4-  -  f(p,T)  -  exp  {a**p  +~-  (=:  -  1 ) }  (19) 

y  o  i  o  1 


Smith  et  al.  [1]  and  Walowit  and  Smith  [2]  used  the  type  I  relation,  and 
they  have  shown  that  for  such  a  relation  a  close  form  formula  for  traction 
computation  in  a  rolling-sliding  contact  may  be  derived.  In  the  case  of 
type  II  relation,  however,  the  integration  is  not  quite  straight-forward  and 
a  numerical  approach  is  preferred,  as  shown  by  Gupta  et  al  [3]. 

The  basic  difference  in  the  two  relations  lies  in  their  physical 
characteristics.  With  the  increased  thermal  heating  during  increasing  slid¬ 
ing  velocity,  the  type  I  relation  gives  a  diminishing  traction  coefficient 
in  the  high  slip  region.  With  the  type  II  relation,  however,  the  thermal 
effect  tends  to  level  off  to  an  asymptotic  value,  which  results  in  a  cor¬ 
responding  steady  value  of  the  traction  coefficient  at  high  slip  velocities. 


Aside  from  the  above  formulation,  some  of  the  more  recent  work  suggests 
alternate  approaches  to  traction  computation  in  a  concentrated  rolling- 
sliding  contact.  Some  of  the  noted  works  are  due  to  Johnson  and  Tevaariverk 
[25],  Bair  and  Winer  [26]  and  Trachman  [6],  The  adaptability  of  these 
models  to  military  lubricants  used  in  rolling  bearing  applications  has  yet 
to  be  investigated. 


2.3  Regression  Analysis  for  Model  Coefficients 

As  discussed  above,  the  determination  of  model  coefficients,  such  as,  u„*, 
a*,  8*,  etc.,  is  an  Important  task  in  the  development  of  a  traction  model. 
In  terms  of  the  above  type  I  or  type  II  relations,  perhaps  the  most  direct 
method  of  determining  the  coefficients  is  to  measure  the  viscosity  as  a 
function  of  both  pressure  and  temperature  and  then  curve-fit  the  data 
to  appropriate  relation.  However,  high-pressure  viscometry  has  been  a  very 
difficult  subject;  therefore,  a  more  serai-empirical  approach,  which  consists 
of  traction  measurement  in  a  single  contact  and  derivation  of  model  coeffi¬ 
cient  from  a  regression  analysis  of  the  traction  data,  has  been  more 
suitable.  Although  such  an  approach  does  not  provide  significant  insight 
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into  the  physics  of  lubricant  behavior,  it  does  result  in  a  realistic  model 
for  computation  of  traction  in  practical  applications. 

The  apparatus  used  for  the  measurement  of  traction  consists  of  two  independ¬ 
ently  driven  disks,  loaded  against  each  other,  when  the  lubricant  is  fed 
through  a  jet.  See  schematic  diagram  in  figure  2.  The  disks  are  normally 
crowned  to  simulate  an  elliptical  contact  of  rolling  bearing  applications. 
The  test  rig  at  the  AFWAL  Aero  Propulsion  Laboratory  (AFWAL/POSL) ,  which  is 
one  of  several  such  facilities,  is  used  in  the  present  investigation. 

Analysis  of  the  experimental  traction  data  consists  of  first  creating 
data  subsets  which  represent  traction-slip  behavior  as  a  function  of  contact 
pressure  under  fixed  rolling  velocity  and  operating  temperature.  Since  a 
reliable  measurement  of  disk  temperature  is  quite  difficult,  the  disk  tem¬ 
perature  is  assumed  to  be  equal  to  the  inlet  temperature.  The  traction 
behavior  is  now  fitted  to  either  the  type  I  or  type  II  relation  and  the 
model  coefficients  are  derived  by  a  least-squares  regression  analysis. 


Consider  the  elliptical  contact  zone  shown  in  figure  2  and  assume  the  slip 
velocity  to  be  constant  over  the  contact  area,  as  will  the  case  in  the  test 
rig.  Using  the  type  I  relation,  equations  (1*1)  to  (16)  may  be  written  as: 


M  -  u o  exp  {a*p  +  0*T © ( 1 -T) }  (20) 

fjj  -  ^^-(l-U)exp  {a*p  ♦  B*T 0 ( 1  — T) )  (21) 

t  -  (  %xp  {a*p  +  B*T o ( 1  -T) )  dU  (22) 

n  JU, 


Straight  forward  integration  of  equation  (21)  and  (22)  shall  result  in  the 
following  expression  for  the  shear  stress: 


where 


A2 


(U2  -  U,)  exp(a*p) 

u„»g»ea<>pU2(U,  -  U,)2 
0K 


slnh  ] A 
/  1  ♦  A2 


(23) 

(2H) 


Since  the  slip  is  constant  over  the  contact  ellipse,  the  total  tractive 
force,  F,  over  the  contact  is  written  as 


F 


f1  f211 

I  abT  d6  rdr 
J0  J0 


2irab 


r 


i  rdr 


(25) 


and  the  global  traction  coefficient. 


K 


F 

Q 


2irab 

Q 


^global’ 

rdr 


written  now  simply  as  k. 


is: 


(26) 


Lubricant 


Figure  2. 


Schematic  Representation  of  a 
Rolling  Disc  Machine  for  Lubricant 
Traction  Measurement 


where  Q  is  the  normal  load  and  a  and  b  are,  respectively,  the  semi-  major 
and  minor  axes  of  the  contaot  ellipse. 

Realizing  that  the  pressure,  p,  in  equations  (23)  and  {2k)  is  elliptical, 
with  a  maximum  Hertz  pressure,  pH, 

p  -  pH  ✓  1  -  r*  (27) 

the  traction  coefficient  may  be  reduced  to  the  following  algebraic  equation: 


<  “  — ; — — 

where 

P  »  a*PH/2.  * 


—^4  fCsinh-1  (Yep)]J 
-  1  [f( fep)  -  f ( V) ] } 

-  U(Ua  -  Ut)  and 


(28) 


f(f)  -  a*. 

Jo  * 

The  integral  f(t)  may  be  computed  numerically,  and  it  may  be  tabulated  over 
the  expected  range  of  Y. 

For  the  purpose  of  regression  analysis,  the  traction  coefficient  is  further 
abbreviated  as 

X  Y  Z  (29) 


3  8Ku„* 

2  PHh  8* 

1 

P  T 

[sinh-1  (fep)]a  -  J-  Cf(fep)  -  f(f)] 

If  the  measured  traction  coefficient  for  any  data  point  i  is  <  .,  then  the 

mi 

squared  deviation  between  the  predicted  and  measured  traction  is  defined  as 


where 


and 


K  - 


X 

Y 


L 


N 


<  )2 
mr 


(30) 


The  three  variables  with  respect  to  which  L  has  to  be  minimized  are  A,  p  and 
V.  It  is  clear  that  for  prescribed  p  and  Y,  the  variation  of  L  with  X  is 
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linear;  therefore,  X  may  be  readily  oomputed  by  conventional  linear  regres¬ 
sion  analysis.  Variations  with  p  and  V  are  nonlinear,  and  for  computing 
these  parameters  the  problem  is  set-up  as: 

f.'s  -  ".!>  sj  ■  v°-  (31) 

and  the  Jacobian 


N  3k,  3k, 

V  {  —L  — i 
-1  1 


♦  (*, 


3*k, 


Kmi^  3x.3x, 


J  k 


(32) 


Once  the  various  derivatives  have  been  expressed  in  terms  of  derivatives  of 
Y  and  Z,  defined  above,  the  set  of  two  simultaneous  equations  (3D  may  be 
solved  by  conventional  Newton-Raphson  techniques. 


A  similar  procedure  can  be  used  for  the  type  II  relation,  but  the  complexity 
of  the  problem  increases  greatly  because  the  entire  problem  is  now  solved 
numerically.  However,  the  coefficients  of  the  type  I  relation  may  be  used 
to  determine  realistic  first  estimates  for  the  type  II  model.  Also,  it  must 
be  noted  that  since  the  pressure-viscosity  coefficient  essentially  in¬ 
fluences  the  tractive  behavior  in  the  low  slip  region,  it  is  quite 
reasonable  to  assume  a*  -  a**.  The  other  two  coefficients  are  computed 
numerically. 


The  model  coefficient,  computed  above,  may  be  tabulated  over  the  range  of 
rolling  speeds  and  operating  temperatures.  Appropriate  values  of  these 
coefficients  may  then  be  used  in  bearing  dynamics  computer  programs,  such 
as  ADORE,  to  determine  the  performance  of  the  bearing. 

In  reality,  the  above  process  is  less  tedious  than  it  appears,  both  in  terms 
of  data  analysis  and  the  subsequent  use  of  the  model  in  bearing  dynamics 
computer  programs.  The  analysis  of  the  data  has  to  be  carried  out  only 
once,  and  all  the  coefficients  may  be  tabulated  for  future  use.  The  use  of 
the  model  in  ADORE  is  quite  straight  forward,  at  least  with  the  type  I 
relation. 


It  must  be  remembered  that  the  computational  efficiency  of  ADORE,  and  other 
bearing  dynamics  programs,  greatly  depends  on  the  complexity  of  the  traction 
model,  since  this  model  operates  in  the  most  inner  loop  of  the  program.  The 
use  of  the  type  II  model,  where  the  shear  stress  is  computed  numerically, 
therefore,  needs  some  caution.  Also,  with  particular  emphasis  on  computa¬ 
tional  simplicity  other  models  [6,25,26]  should  be  further  explored.  The 
regression  analysis  described  above  provides  a  generic  approach  for  the 
computation  of  coefficient  of  any  arbitrary  model  from  the  experimental 
traction  data  obtained  over  the  expected  range  of  operating  conditions. 


3.  RESULTS 


Based  on  the  analytical  formulation  described  in  the  preceding  section, 
experimental  traction  data  forms  the  basic  foundation  to  the  traction  model 
development.  For  the  present  investigation,  the  data  obtained  at  the  AFWAL 
Aero  Propulsion  Laboratory  for  the  MIL-L-7808  and  MIL-L-23699-type  military 
lubricants  is  considered,  and  the  various  coefficients  for  the  type  I  and 
type  II  viscosity-pressuhe-temperature  relations  are  derived.  The  model  is 
then  incorporated  in  the  bearing  dynamics  computer  program,  ADORE,  and  the 
general  slip  behavior  of  a  roller  is  simulated  in  a  cylindrical  roller 
bearing. 

Most  of  the  available  data  were  obtained  with  a  set  of  76.20  mm  (3.0  inch) 
diameter  disks  with  a  crown  radius  of  91  M.40  mm  (36  inch);  in  one  of  the 
data  sets  the  crown  radius  was  set  equal  to  the  rolling  radius  to  simulate  a 
circular  point  contact.  The  inlet  oil  temperature  was  varied  in  the  range 
of  290  to  366  °K;  the  rolling  velocity  was  varied  from  21  to  M5  M/Sec;  and 
the  contact  pressures  ranged  from  5M0  to  860  MPa  with  the  91M.M0  mm  crowned 
disks  and  1 MM0  to  2050  MPa  in  the  case  of  a  circular  point  contact.  In  all 
cases  the  lubricant  was  fed  through  a  Jet  and  the  inlet  temperature  was 
closely  monitored.  In  3ome  cases  measurement  of  disk  surface  temperature 
was  also  attempted. 


3.1  General  Traction  Behavior 

The  general  traction  behavior  of  most  lubricants  in  an  elastohydrodynamic 
contact  is  described  in  figure  3a.  The  traction  coefficient  first  in¬ 
creases,  almost  linearly,  with  slip,  it  peaks  to  a  maximum  value,  and  then 
with  the  increasing  thermal  effect  with  higher  slip  velocities  it  begins  to 
drop.  The  discussion  presented  in  the  preceding  section  clearly 
demonstrates  that  the  analytical  model  considered  herein  preserves  such  a 
general  behavior.  In  fact,  figure  3a  shows  the  fit  of  the  model  to  actual 
traction  data  obtained  at  certain  experimental  conditions.  See  Appendix  A 
for  the  explanation  of  the  various  symbols  used  in  the  computer-generated 
traction  plots,  such  as,  figure  3a,  and  others,  throughout  this  report.  The 
actual  experimental  points  are  marked  with  a  square,  while  the  solid  line 
represents  prediction  of  the  model  based  on  the  type  I  relation.  The 
results  with  the  type  II  relation  are  similar.  See  figure  3b.  Clearly,  the 
fit  between  model  predictions  and  experimental  data  is  very  good.  Thus  for 
any  given  traction-slip  curve,  the  model  coefficients  may  be  derived  to 
provide  a  close  simulation  of  experimental!/  observed  traction  behavior. 

Since  the  practical  significance  of  a  traction  model  lies  in  its  predictive 
capabilities  over  a  reasonably  wide  range  of  operating  conditions,  the  close 
fit  of  the  data  shown  in  figures  3a  and  3b,  although  necessary,  does  not 
provide  sufficient  confidence  in  the  model  predictions.  Clearly  the  deriva¬ 
tion  and  tabulation  of  model  coefficients  for  each  operating  condition  are 
impractical.  Realistic  objective  of  a  model  is  to  simulate  the  observed 
traction  behavior  over  the  entire  range  of  data  with  one  set  of  model 
coefficients.  With  such  an  ultimate  objective,  the  regression  analysis  is 
carried  out  over  the  range  of  contact  pressures  at  prescribed  rolling 
velocity  and  inlet  temperature.  Although,  the  coefficients  so  derived  may 
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Figure  3a.  Typical  traction  behavior  and  the  predictions  of  the  type  I  model 


Figure  3b.  General  traction  behavior  and  the  predictions  of  the  type  II  model 


be  further  curve-fitted  to  cover  the  variation  over  rolling  velocity  and 
temperature,  they  have  a  substantial  significance  in  rolling  bearings  where, 
for  a  given  application,  the  variations  in  rolling  velocities  and  operating 
temperatures  is  not  very  large.  Thus,  as  a  first  step  to  traction  model 
development,  the  model  coefficients  in  the  present  investigation  are  tabu¬ 
lated  over  the  rolling  velocities  and  inlet  oil  temperatures. 


3.2  Traction  Sensitivity  to  Model  Coefficients 

In  order  to  establish  some  physical  significance  of  each  of  the  model  coef¬ 
ficient,  it  is  necessary  to  investigate  the  sensitivity  of  traction  behavior 
to  each  model  coefficient  before  presenting  the  results  of  the  regression 
analysis  carried  out  in  the  present  investigation. 

There  are  essentially  three  parameters  of  the  models  considered  herein:  the 
reference  viscosity,  y0*  or  y0**,  the  pressure  viscosity  coefficient,  a*  or 
a**;  and  the  temperature-viscosity  coefficient,  B*  or  B**.  Figures  4a  and 
4b  show  the  variation  in  traction  as  a  function  of  the  reference  viscosity 
for  the  type  I  and  type  II  models  respectively.  As  expected,  the  increase 
in  viscosity  results  in  a  corresponding  increase  in  traction  for  both  types 
of  models.  The  results  of  the  pressure-viscosity  coefficients  are  also 
similar,  as  seen  in  figures  5a  and  5b.  The  temperature-viscosity  parameter 
demonstrates,  perhaps,  the  most  interesting  results;  traction  as  a  function 
of  B*  is  shown  in  figure  6a  while  similar  results  as  a  function  of  B**  for 
the  type  II  model  are  shown  in  figure  6b.  Note  that  at  high  slip  velocities 
the  type  I  model  results  in  a  relatively  faster  drop-off  of  traction  coeffi¬ 
cient  compared  to  the  type  II  model;  this  is  primarily  due  to  the  thermal 
behavior  of  the  two  constitutive  relations.  At  low  slip,  however,  the  two 
models  are  essentially  identical  when  the  pressure-viscosity  coefficients 
for  the  two  models  are  the  same.  In  fact,  it  is  the  pressure-viscosity 
coefficient  which  essentially  establishes  the  slope  of  the  traction  curve  at 
low  slip  velocities,  as  discussed  earlier. 

The  overall  difference  between  the  type  I  and  type  II  relations  is,  there¬ 
fore,  only  significant  at  high  slip  velocities.  Under  excessive  slip,  as 
the  temperature  of  the  lubricant  increases,  the  effective  viscosity  con¬ 
tinues  to  drop  at  a  steady  rate  for  the  type  I  model;  therefore  the  traction 
coefficient  continues  to  drop.  In  the  type  II  model,  on  the  other  hand,  the 
rate  at  which  the  viscosity  drops  gradually  slows  down  due  to  the  inverse 
temperature  effect;  this  results  in  a  somewhat  asymptotic  behavior  of  the 
traction  coefficient  at  very  high  slip  velocities.  Such  a  difference  in 
behavior  is  seen  in  figure  7,  where  traction  predictions  by  the  two  models 
are  shown  under  Identical  operating  conditions. 


3.3  MIL-L-23699-Type  Lubricant 

The  ambient  properties  of  the  MIL-L-23699-type  lubricant  are  well  documented 
in  the  literature.  In  the  present  investigation  the  data  compiled  by  Gupta 
et  al.  [3]  is  used  for  the  computation  of  film  thickness.  Also,  based  on 
the  model  correlations  reported  in  this  work  the  pressure-viscosity  coeffi¬ 
cient  in  the  high  pressure  contact  region,  a*  or  a#*,  is  assumed  to  be  5.8E- 
09  1/Pa. 
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Figure  4a.  The  influence  of  reference  viscosity,  /u*  (MU),  in  the 
type  I  model. 
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Figure  4b.  The  influence  of  reference  viscosity,  u**  (MU),  in  the 
type  II  model. 
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Figure  5a.  The  influence  of  pressure-viscosity  coefficient,  a  (ALFA), 


Figure  5b.  The  influence  of  pressure- viscosity  coefficient,  a  (ALFA), 
in  the  type  II  model. 
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Figure  6a.  The  influence  of  temperature-viscosity  coefficient,  0  (BETA), 
in  the  type  l  model. 


Figure  6b  The  influence  of  temperature- viscosity  coefficient,  |8**(BETA), 
in  the  type  II  model. 


Figure 


Figures  8a  and  8b  show  the  experimental  traction  data  (marked  points)  and 
model  predictions  (solid  lines)  for  the  type  I  and  type  II  models  respec¬ 
tively,  at  an  inlet  temperature  of  31*1*  °K  and  at  a  rolling  speed  of  30.5 
M/Sec.  Clearly,  the  fit  between  the  data  and  model  predictions  is  very  good 
for  both  types  of  viscosity-pressure-temperature  relations.  The  fit  is  also 
very  good  at  the  higher  rolling  speeds  and  higher  inlet  temperatures, 
figures  9a  and  9b  show  the  results  at  366  °K  inlet  temperature  and  at  a 
rolling  speed  of  H5.5  M/Sec,  which  represents  the  extreme  operating  condi¬ 
tions  of  the  available  data.  At  the  low  end  of  inlet  temperatures  and 
rolling  speeds,  however,  there  is  substantial  difference  between  model 
predictions  and  experimental  observations.  As  shown  in  figure  10a,  the 
regression  coefficients  for  the  type  I  model  provide  a  relatively  good  fit 
to  the  medium  pressure  data,  it  under-predicts  traction  at  higher  pressures 
and  the  predictions  are  higher  than  the  experimental  observations  at  lower 
pressures.  With  the  type  II  relation,  however,  the  regression  analysis 
results  in  a  fairly  good  fit  at  higher  pressures,  while  the  predictions  are 
higher  than  the  experimental  observations  at  lower  pressures,  as  shown  in 
figure  10b. 

On  the  whole  correlations  between  the  model  predictions  and  experimental 
observations  are  equally  good  for  both  the  type  I  and  type  II  models.  Plots 
for  the  entire  data  set  are  included  with  the  model  predictions  in  Appendix 
A  and  the  model  coefficients,  computed  from  the  regression  analysis,  are 
summarized  in  table  1 . 

A  close  examination  of  the  coefficients  tabulated  in  table  1  reveals  that 
the  reference  viscosity  parameters,  u0*  and  u0**»  generally  decrease  with 
both  rolling  speed  and  temperature.  This  is  in  complete  agreement  with  the 
trends  observed  earlier  by  Gupta  et  al  [3].  The  variation  in  8*,  however, 
unlike  the  earlier  data  [3],  does  not  appear  to  be  significant;  this  sug¬ 
gests  a  further  refinement  in  the  model  by  carrying  out  the  regression 
analysis  over  the  viscosity  parameter  only.  For  the  type  II  model,  however, 
there  is  a  noticeable  variation  in  8**  as  a  function  of  rolling  speed  and 
temperature;  it  seems  to  decrease  with  rolling  speed  and  increase  with 
temperature.  Such  a  trend  is  in  agreement  with  the  earlier  data  [3].  A 
rigorous  quantitative  comparison  of  the  coefficients  derived  by  Gupta  et  al 
[3]  to  those  obtained  under  the  present  investigation  is  not  quite 
straightforward  because  the  disk  temperature  data  is  not  available  in  the 
current  data  base.  However,  in  view  of  the  present  observation  that,  at 
least  for  the  type  I  relation,  the  traction  model  may  be  reduced  to  only  a 
single  viscosity  parameter,  it  will  be  interesting  to  repeat  the  regression 
analysis  of  both  the  earlier  and  the  current  data  with  a  fixed  8*  and  with 
the  disk  temperatures  equated  to  the  inlet  temperature. 


3.^  MIL-L-7808-Type  Lubricant 

The  MIL-L-7808-type  lubricant  is,  perhaps,  the  most  widely  used  military 
lubricant  and  its  rheological  behavior  has  been  studied  by  a  number  of 
investigators.  In  the  present  investigation,  the  ambient  properties  docu¬ 
mented  by  Gupta  et  al.  [3]  are  used  for  the  film  thickness  computation  and 
also,  based  on  the  parametric  study  reported  in  this  work  the  pressure- 
viscosity  coefficient  in  the  high-pressure  contact  region  Is  assumed  to  be 
5.22E-09  1/Pa,  as  a  baseline. 
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Figure  8a.  Fit  of  the  type  I  model  to  the  MIL-L-23699  data  at  344  K  and 
at  a  rolling  speed  of  30.5  meters/second. 


Figure  9b.  Fit  of  the  type  II  nodel  to  the  MIL-L-23699  data  at  366  K  and  at  rolling  speed  of 
45.5  meters/second . 
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Figure  10b.  Fit  of  the  type  I  model  to  the  MIL-L-23699  data  at  294  °K  and  at  a  rolling  speed 
of  21.15  meters/second. 


Table  1 


Model  Coefficients  for  the  MIL-L-23699  Lubricant 
Cor  elations  at  Inlet  Temperature 


a*  -  a** 

-  5.802  x10~9 

1/Pa 

Temp 

Roll  Velo 

8* 

( °K) 

(M/S) 

(Pa.S) 

(1/#K) 

(Pa.S) 

(°K) 

294 

21.146 

0.2995 

0.04945 

0.4842 

5611 

294 

30.482 

0.2651 

0.05075 

0.2771 

5110 

294 

45.544 

0.07182 

0.03998 

0.09958 

4348 

3  22 

21.146 

0.1290 

0.05621 

0.1571 

6851 

322 

30.482 

322 

45.544 

0.04668 

0.04984 

0.05358 

5491 

344 

21.146 

0.04634 

0.04905 

0.05646 

7632 

344 

30.482 

0.03917 

0.04378 

0.04319 

6079 

344 

45.544 

0.03103 

0.04692 

0.03215 

5969 

366 

21.146 

0.01915 

0.06207 

0.02058 

9861 

366 

30.482 

0.01752 

0.04668 

0.01896 

7307 

366 

45.544 

0.01625 

0.04643 

0.01498 

6113 

The  model  correlations  are  very  similar  to  those  observed  for  the  MIL-L- 
23699  type  lubricant.  The  fit  between  the  data  and  model  predictions  is 
extremely  good  throughout  the  entire  range  of  operating  conditions,  except 
for  the  few  data  sets  obtained  at  the  lowest  inlet  temperature.  Figures  11a 
and  11b  show  the  fits  at  an  inlet  temperature  of  322  °K  and  at  the  high  end 
of  rolling  speeds,  45.5  M/Sec,  for  the  type  I  and  type  II  models 
respectively.  The  fit  at  the  higher  temperature  of  344  °K  is  shown  in 
figures  12a  and  12b,  respectively  for  the  two  types  of  viscosity-pressure- 
temperature  relations.  At  the  lower  inlet  temperature,  294  °K,  similar  to 
the  MIL-L-23699  type  lubricant,  the  type  I  model  seems  to  better  fit  the 
data  at  medium  pressures,  while  the  type  II  relation  fits  better  at  higher 
pressures,  as  seen  in  figures  13a  and  13b.  A  worthy  observation  here  is  the 
fact  that  there  is  a  definite  bias  in  this  set  of  experimental  data,  i.e., 
the  traction  coefficients  at  the  positive  slip  velocities  are  different  from 
those  at  the  corresponding  negative  value  of  slip  velocities,  particularly 
at  the  high  end  of  contact  pressures.  One  of  the  factors  contributing  to 
this  difference  in  traction  is  the  different  rolling  speed  at  the  cor¬ 
responding  points  on  the  positive  and  negative  sides  of  the  data,  due  to  the 
fact  that  the  variation  in  slip  is  obtained  by  varying  the  velocity  of  only 
one  of  the  disks.  However,  the  difference  seen  in  figures  13a  and  13b, 
cannot  be  fully  accounted  for  by  the  varying  rolling  speed  alone.  Perhaps, 
the  torque  measurement  and  calibration  systems,  and  the  general  nature  of 
power  dissipation  in  the  support  bearings,  particularly  at  lower  oil  tem¬ 
peratures,  need  further  examination. 

Based  on  the  above  discussion  it  may  be  concluded  that  except  for  the  data 
obtained  at  the  low  end  of  inlet  temperatures,  there  is  a  good  agreement 
between  the  model  predictions  and  experimental  observations.  The  computer 
plots  for  the  entire  range  of  data  along  with  model  predictions,  for  both 
the  type  I  and  type  II  models,  are  included  in  Appendix  A.  The  model  coef¬ 
ficients,  computed  from  the  regression  analysis,  are  summarized  in  table  2. 

The  trends  of  the  various  model  parameters  are  similar  to  those  discussed 
above  for  the  MIL-L-23699  lubricant  and  they  are  in  good  agreement  to  those 
obtained  earlier  by  Gupta  et  al.  [3],  except  that  the  variation  in  the 
temperature-viscosity  parameter,  0#,  is  not  significant  in  the  current  datr. 
Based  on  this  finding  it  is  necessary  to  repeat  the  regression  analysis  with 
fixed  0*  for  both  the  earlier  and  current  data.  Also,  the  disk  temperature 
in  the  earlier  data  should  be  set  equal  to  the  inlet  temperature  for  a  one 
to  one  comparison  with  the  current  data. 


3.5  Effect  of  Contact  Geometry 

All  of  the  experimental  data  discussed  above  was  obtained  with  u  set  of 
crowned  disk  specimens,  with  a  rolling  radius  of  38.10  mm  and  a  crown  radius 
of  914.40  mm.  Thus  the  shape  of  the  contact  zone  was  a  relatively  narrow 
ellipse  which  simulates  either  a  closely  conformed  ball/race  contact  in  a 
ball  bearing  or  a  roller/race  line  contact  in  cylindrical  roller  bearings. 
With  the  objective  of  investigating  the  influence  of  contact  geometry  one 
set  of  traction  data,  for  the  MIL-L-7808-type  lubricant,  was  obtained  with 
spherical  disks,  where  the  crown  radius  is  set  equal  to  the  rolling  radius; 
thus,  a  circular  point  contact  is  simulated. 
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Figure  12a.  Fit  of  the  type  I  model  to  the  MIL-L-7808  data  at  344  K  and  at  a  rolling  speed  of 
45.5  meters/second. 
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Table  2 


Model  Coefficients  for  the  MIL-L-7808  Lubricant 


Correlations  at  Inlet  Temperature 


a*  -  a** 

-  5.221  xio"9 

1/Pa 

Temp 

Roll  Velo 

Mo* 

6* 

Mo** 

0«» 

( °K) 

(M/S) 

(Pa.S) 

( 1/*K) 

(Pa.S) 

( *K) 

294 

21.1146 

0.371414 

0.014397 

0.6313 

*4932 

29** 

30.H82 

0.2289 

0.04243 

0.14060 

14812 

294 

145.51414 

0.1013 

0.03855 

0.1558 

14187 

322 

21.1H6 

0.10814 

0.0143141 

0.1591 

6117 

322 

30.1482 

0.08075 

0.03579 

0.1047 

H7143 

322 

145.51414 

0.07136 

0.03138 

0.071176 

3818 

344 

21.H46 

0.014119 

0.05071 

0.014771 

7231 

344 

30.1482 

0.03073 

0.035514 

0.03712 

51448 

31414 

145.514H 

0.02762 

0.03130 

0.03032 

4405 

Figures  14a  and  1 4b  show  a  typical  set  of  traction  data  and  its  correlations 
with  the  type  I  and  type  II  models.  It  Is  seen  that  the  type  I  model  fits 
the  experimental  behavior  fairly  well,  although  the  fit  is  not  as  good  as 
that  observed  earlier  for  the  elliptical  contact  data.  Note  the  rather  high 
slope  of  the  traction  curve,  which  results  from  the  relatively  high  contact 
stresses  in  the  point  contact  data.  Under  the  very  high  contact  stress  the 
predictions  of  the  type  II  model  are  somewhat  questionable,  particularly 
with  regard  to  computation  of  thermal  reduction  factor  for  the  lubricant 
film,  the  viscous  shear  heating  in  the  high-pressure  contact  zone,  and  the 
positive  traction  slope  in  the  high  slip  region,  as  seen  in  figure  14b. 

The  computed  model  coefficients  for  this  set  of  data  are  tabulated  in  table 
3  and  the  traction  curves  are  included  in  Appendix  A.  Similar  to  the  ellip¬ 
tical  contact  data,  discussed  above,  the  coefficients  in  table  3  are 
computed  with  the  assumption  that  the  surface  temperature  of  the  disks  is 
equal  to  the  inlet  temperature,  The  general  trends  of  the  coefficients  are 
quite  similar  to  those  discussed  above  for  the  line  contact  data.  However, 
in  view  of  the  substantially  increased  scatter  between  the  experimental  data 
and  model  predictions,  particularly  in  terms  of  the  very  high  traction 
slopes  at  low  slip  rates  and  very  high  contact  pressures,  it  is  necessary  to 
modify  the  model  to  allow  for  viscoelastic  effects,  after  ensuring  the 
repeatability  of  the  data.  Such  refinements  may  be  based  on  the  more  recent 
work  by  Johnson  and  Tevaariverk  [26],  and  by  Bair  and  Winer  [27]. 

Since  the  point  contact  data  includes  some  measurement  of  disk  tempera¬ 
tures,  the  regression  analysis  is  repeated  with  the  actual  disk 
temperatures;  the  resulting  coefficients  are  tabulated  in  table  4.  Although 
the  scatter  between  the  experimental  data  and  the  model  predictions,  see  the 
data  set  plots  in  Appendix  A,  is  relatively  unchanged,  there  is  a  rather 
erratic  change  in  the  model  coefficients.  The  jump  in  the  viscosity 
parameters  at  322  °K,  in  table  4,  suggests  possible  experimental  dis¬ 
crepancies  in  the  temperature  data. 


3.6  Rolling  Bearing  Performance  Modeling 

The  significance  of  the  traction  model  lies  In  the  prediction  of  component 
behavior  and  the  simulation  of  practical  mechanical  systems  where  the 
lubricant  is  used  in  practice.  For  the  purpose  of  illustration,  the  model 
coefficients  computed  above  for  the  type  I  model  are  input  to  the  bearing 
dynamics  computer  model,  ADORE,  and  the  slip  characteristics  of  a  roller  In 
a  cylindrical  roller  bearing  are  investigated  as  the  roller  travels  through 
the  load  zone  in  the  radial  loaded  bearing.  Figure  15  shows  the  length  of 
simulation  and  the  normal  load  variation  as  the  roller  travels  from  a  rela¬ 
tively  unloaded  state  to  the  maximum  load  state  and  then  exits  out  of  the 
load  zone. 

As  the  roller  travels  through  such  varying  contact  loads,  the  tractive 
forces  at  the  roller/race  interaction,  which  are  derived  from  the  lubricant 
traction  model,  determine  the  acceleration  of  the  roller;  hence  roller  slip 
or  skid  is  very  closely  coupled  to  the  lubricant  traction  behavior.  Figure 
16a  shows  the  simulated  behavior  under  the  type  I  model  at  294  °K.  The 
operating  conditions  at  the  roller/race  interface,  particularly  the  rolling 
speed,  are  within  the  range  of  experimental  data  used  to  derive  the  model 
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Typical  correlation  of  the  type  I  model  to  the  high-pressure  circular  contact 
traction  data  obtained  with  the  MIL-L-7808  lubricant. 
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Figure  14b.  Typical  correlation  of  the  type  II  model  to  the  high-pressure  circular  contact 
traction  data  obtained  with  the  MIL-L-7808  lubricant. 


Table  3 


Model  Coefficients  for  the  H!L-L-7808  Lubricant 
Point  Contact  Data 
Correlations  at  Inlet  Temperature 


.  a#* 

-9 

-  5.221  xIO 

1/Pa 

Temp 

(°K) 

Roll  Velo 
(M/S) 

Wo* 

(Pa.S) 

8* 

(1/°K) 

Wo** 

(Pa.S) 

e«* 

( °K) 

294 

294 

294 

21.146 

30.482 

45.544 

0.1317 

0.03520 

0.02477 

0.04559 

0.01755 

0.02878 

0.06474 

0.01925 

0.005355 

4019 

2273 

2832 

322 

322 

322 

21 .146 
30.482 
45.544 

0.07780 

0.06522 

0.06168 

0.05270 

0.04312 

0.04114 

0.03988 

0.02796 

0.01542 

5305 

4132 

3420 

344 

21.146 

344 

344 

30.482 

45.544 

0.01655 

0.03636 

0.009870 

4611 

Table  4 

Model  coefficients  for  the  MIL-L-7808  Lubricant 


Point  Contact  Data 
Correlations  at  Disk  Temperature 


m  a** 

-  5.802  xIO 

1/Pa 

Temp 

(°K) 

Roll  Velo 
(M/S) 

Wo* 

(Pa.S) 
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Figure  15,  Typical  roller/race  load  variation  in  a  radially  loaded 
cylindrical  roller  bearing. 
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coefficients.  As  the  load  on  the  roller  increases,  it  demands  a  definite 
traction  force  from  the  lubricant  in  order  to  maintain  a  rolling  configura¬ 
tion  at  both  races.  Thus  the  slip  on  the  driving  race,  which  is  the  inner 
race  in  the  present  case,  begins  to  increase  in  search  of  increasing 
traction;  the  slip  will  normally  increase  until  the  required  acceleration 
is  met.  With  a  liquid-lubricated  bearing,  the  increasing  slip  results  in  a 
reduction  in  film  thickness  due  to  shear  heating  of  the  lubricant.  This  in 
turn  results  in  a  reduction  in  tractioi.  through  the  film,  but  it  increases 
the  probability  of  metal  contact  as  the  film  thickness  approaches  the  rms 
surface  roughness  of  the  interacting  surfaces.  Such  a  mechanism  for  the 
MIL-L-7808  lubricant  is  clearly  seen  in  figure  16a.  It  is  seen  that  since 
the  lubricant  traction  force  is  always  less  than  that  required  to  keep  the 
roller  in  a  rolling  state,  the  slip  at  the  roller  to  inner  race  interface 
increases  to  a  maximum  when  the  roller  reaches  the  maximum  load.  The  ther¬ 
mal  effect  in  the  lubricant  film  through  this  slippage  reduces  the  film 
thickness  to  almost  zero  under  the  maximum  load.  Note  that  since  the  solu¬ 
tions  shown  in  figure  16a  do  not  include  any  metal  contact  and  are  strictly 
based  on  fluid  traction,  the  contact  losses,  particularly  at  the  inner  race, 
may  be  greatly  under  estimated. 

If  the  operating  temperature  is  increased  to  3^  °K,  the  upper  limit  of 
experimental  data,  both  the  lubricant  film  thickness  a-  d  traction  further 
reduce  but  the  general  behavior  of  roller  remains  unchanged.  This  is  shown 
in  figure  16b,  which  is  qualitatively  similar  to  figure  16a. 

A  practical  question  which  arises  from  the  above  discussion  concerns  the 
determination  of  optimum  lubricant  behavior  for  acceptable  bearing 
performance.  It  is  necessary  that  the  lubricant  provide  adequate  traction 
to  prevent  excessive  skid  on  the  roller,  yet  the  traction  should  be  small 
enough  to  keep  the  heat  generation  in  the  bearing  within  reasonable  limits. 
A  realistic  modeling  of  lubricant  behavior  is,  therefore,  extremely 
important.  To  establish  this  point  a  bit  further,  consider  the  traction 
curve  of  figure  17,  which  provides  somewhat  larger  tractions  than  those 
observed  with  the  MIL-L-7808  lubricant.  Under  such  a  behavior,  the  traction 
forces,  as  the  roller  begins  to  slip,  rapidly  meet  the  roller  acceleration 
requirements;  and  since  the  traction  response  at  the  outer  and  inner  races 
is  identical  under  the  hypothetical  model,  the  roller  slip  rates  and  trac¬ 
tion  coefficients  adjust  such  that  the  heat  generated  at  the  outer  and  inner 
race  interactions  are  almost  equal.  Such  a  behavior  is  shown  in  figure  18. 
It  is  also  seen  that  the  slip  rates  in  figure  18  are  lower  than  those  shown 
in  figures  16a  and  16b;  however,  the  heat  generated  at  roller/race  contact 
is  significantly  larger  due  to  higher  traction  coefficients,  which  implies 
that  certain  optimization  of  the  lubricant  properties  is  necessary.  Such  an 
optimization  can  be  easily  performed  by  carrying  out  a  parametric  evaluation 
of  the  roller  motion  as  a  function  of  the  coefficients  of  the  traction 
model.  Thus  the  development  of  a  traction  model  is  significant  for  the 
design  and  performance  simulation  of  mechanical  components,  such  as  rolling 
bearings. 
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Figure  17.  A  hypothetical  traction-slip  curve. 
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Figure  18.  The  simulation  of  roller/race  slip  in  a  cylindrical  roller  bearing 
with  the  hypothetical  traction  model  shown  in  figure  17. 
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4.  CONCLUSIONS 


The  feasibility  of  a  traction  model,  based  on  viscosity-pressure-temperature 
oonstitutive  relations  in  the  high-pressure  zone  of  elastohydrodynamic  con¬ 
tacts,  is  demonstrated  for  two  military  lubricants  meeting  MIL-L-23699  and 
MIL-L-7808  specifications.  The  model  predictions  are  shown  to  agree  fairly 
well  with  the  experimental  data  obtained  at  the  AFWAL  Aero  Propulsion 
Laboratory,  with  cylindrical  disk  specimens  resulting  in  an  elliptical 
contact. 

Among  the  two  types  of  viscosity-pressure-temperature  relations,  type  I  with 
exponential  variation  in  viscosity  with  pressure  and  temperature,  and  type 
II  with  the  viscosity  varying  exponentially  with  pressure  and  inverse  tem¬ 
perature,  the  regression  coefficients  for  the  type  I  model  indicate  that  the 
entire  variation  may  be  represented  by  a  variation  in  a  reference  viscosity 
parameter;  thus  a  greatly  simplified  relation  for  traction  may  be  obtained. 
Also,  since  the  mathematical  simplicity  of  the  type  I  relation  allows  for  an 
effective  incorporation  of  the  traction  model  in  sophisticated  component 
models,  such  as  the  advanced  bearing  dynamics  computer  code  ADORE,  it  offers 
significant  potential  for  the  performance  simulation  of  practical  components 
as  a  function  of  lubricant  behavior  and,  therefore,  the  development  of 
lubricants  for  advanced  applications. 

In  the  case  of  very  high  pressure  circular  contacts,  noticeable  scatter 
between  the  experimental  data  and  model  predictions  is  observed,  par¬ 
ticularly  in  terms  of  modeling  the  rather  high  traction  slope  in  the  very 
low  slip  region.  This,  to  some  extent,  suggests  a  viscoelastic  behavior  of 
the  lubricant.  Further  model  refinements  to  allow  for  such  physical  be¬ 
havior  are  therefore  necessary  for  applications  with  extremely  high  contact 
pressures. 

The  close  coupling  between  the  lubricant  behavior  and  component  performance 
is  demonstrated  by  modeling  the  roller  slip  as  a  function  of  lubricant  be¬ 
havior  in  a  cylindrical  roller  bearing.  It  is  demonstrated  that  the 
development  of  a  realistic  traction  model  is  vital  to  the  selection  of  op¬ 
timum  lubricant  behavior  for  ?  given  application. 


5.  RECOMMENDATIONS  FOR  FUTURE  DEVELOPMENTS 


The  findings  of  the  present  investigation  demonstrate  the  feasibility  of  the 
overall  approach  to  traction  model  development  through  a  least-squared 
regression  fit  of  a  physical  model  to  actual  traction  data.  However,  the 
ultimate  development  of  a  model  for  military  applications  requires  sig¬ 
nificant  further  development.  Based  on  the  conclusions  of  the  present  work, 
some  of  the  specific  tasks  requiring  further  development  are  listed  below: 

1.  The  regression  analysis  for  the  type  I  relation  should  be  repeated  for 
the  present  data  with  the  entire  variation  in  traction  modeled  in  terms 
of  the  variation  of  the  reference  viscosity  parameter.  A  similar 
analysis  should  be  carried  out  for  other  sets  of  data  available  for 
similar  lubricants.  If  such  an  analysis  results  in  acceptable 
correlations  then  a  significant  simplification  in  the  rheological  model 
may  be  accomplished,  and  this  will  contribute  to  a  substantial 
improvement  in  the  computing  effort  required  to  model  component  behavior 
as  a  function  of  lubricant  properties. 

2.  The  traction  model  should  be  extended  to  include  viscoelastic  effects. 
Appropriate  modifications  to  the  regression  analysis  should  then  be  made 
to  allow  for  the  computation  of  model  coefficients  from  laboratory 
traction  data.  This  will  significantly  improve  traction  predictions  at 
very  high  contact  pressures  and  in  the  low  slip  region. 

3.  Although  the  Isothermal  film  thickness  can  be  easily  computed  from  the 
available  formulae,  the  computation  of  thermal  reduction  factors  becomes 
questionable  as  the  available  numerical  solutions  are  extrapolated  to 
very  high  slip  velocities,  which  are  sometimes  encountered  in  roller 
bearings.  Also,  in  ball  bearings,  when  the  ball/cage  collisions  are 
severe,  excessive  ball/race  slip  or  skid  is  encountered.  Thus  some 
refinements  in  the  thermal  model  for  the  elastohydrodynamic  film 
thickness  are  necessary  to  Improve  the  predictions  at  high  slip 
velocities. 

4.  Some  of  the  very  recent  research  has  established  that  the  lubricant  film 
thickness  in  elastohydrodynamic  contacts  should  also  be  modified  to 
allow  for  the  surface  roughness  of  the  interacting  surfaces.  The 
development  of  "roughness  factors",  based  on  the  available  numerical 
solutions,  is  therefore,  proposed. 

5.  Since  both  the  film  thickness  and  lubricant  behavior  in  the  high- 
pressure  contact  zone  play  an  important  role  in  the  prediction  of 
traction,  it  will  be  desirable  to  carryout  simple  experiments  where  the 
film  thickness  and  traction  may  be  measured  3imalteneously.  The  optical 
rig  at  the  AFWAL  Aero  Propulsion  Laboratory  may  be  used  to  obtain  such 

experimental  data. 

6.  Fundamental  measurements  of  lubricant  rheology  are  vital  to  realistic 
physical  interpretation  of  the  model  coefficients.  For  example, 
significant  physical  Justification  to  the  pressure-viscosity 
coefficients,  derived  from  regression  analysis  of  the  traction  data,  may 
be  obtained  if  the  high-pressure  viscosity  data  is  available.  Indeed, 
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high-pressure  viscometry  is  a  difficult  task,  however,  the  recent 
developments  at  the  AFWAL  Aero  Propulsion  Laboratory  provide 
substantial  advancements  of  such  techniques.  It  is,  therefore, 
suggested  that  the  model  coefficients  should  be  compared  to  the  high- 
pressure  viscoraetry  data  and  appropriate  correlations  should  be 
developed. 

7.  The  maximum  rolling  speeds  and  operating  temperatures  at  which  most  of 
the  available  traction  data  was  obtained  are  well  below  the  actual 
conditions  in  high-speed  rolling  bearings  used  in  a  wide  range  of 
military  applications.  The  models  developed  on  the  basis  of  such  data 
have  to  be  greatly  extrapolated  when  they  are  actually  used  in  bearing 
performance  modeling.  Experimental  work  to  extend  the  operating  range 
of  the  traction  experiments  is  therefore  necessary.  Perhaps  traction 
measurements  with  larger  diameter  disks,  compared  to  the  current  76  mm 
diameter  disks,  will  be  desirable  if  the  maximum  disk  speed  is  an 
experimental  limitation. 

8.  The  recent  traction  data  obtained  at  the  AFWAL  Aero  Propulsion 
Laboratory  for  heavily  loaded  circular  point  contact  Is,  perhaps,  the 
first  data  set  of  its  kind.  Such  data  are  indeed  necessary  over  a  wider 
range  of  contact  geometry  in  order  to  establish  the  effect  of  the 
ellipticity  parameters  on  traction  and  resulting  behavior  of  a  bearing. 
Also,  some  of  traction  experiments  with  varying  geometry  may  be 
carried  out  under  controlled  misalignment  to  simulate  the  effect  of  spin 
velocities  normal  to  the  contact.  The  present  traction  rig  at  the  AFWAL 
Aero  Propulsion  Laboratory  apparently  has  such  a  capability;  therefore, 
such  experiments  may  not  be  very  difficult. 

9.  After  establishing  adequate  validation  of  a  traction  model  on  the  basis 
of  actual  traction  data  obtained  with  single  contacts,  it  is  essential 
to  prove  the  significance  of  the  model  in  component  modeling.  The 
modeling  of  bearing  performance  for  prescribed  traction  models,  and 
comparison  of  the  bearing  performance  predictions  to  bearing  test  data, 
is  therefore  necessary.  Such  test  data  has  been  obtained  for  high-speed 
ball  bearings  at  the  AFWAL  Aero  Propulsion  Laboratory.  However,  due 
to  a  combined  slip  and  spin  of  the  balls,  the  kinematics  of  a  ball/race 
contact  is  extremely  complex,  and  a  validation  of  the  traction  model 
for  such  bearings  requires  extensive  analytical  modeling  of  the  contact. 
The  validation,  therefore,  results  in  a  number  of  questions  with  regard 
to  both  contact  modeling  and  applicability  of  traction  models  developed 
from  greatly  simplified  contacts.  Perhaps  experiments  with 
cylindrical  roller  bearings,  where  the  roller/race  contacts  are 
practically  free  of  any  spin  normal  to  the  contact  and  the  slip  is 
fairly  uniform  over  the  contact,  at  least  for  carefully  aligned 
bearings,  may  provide  validation  of  a  traction  model  on  a  more 
fundamental  level.  Once  such  a  validation  of  the  model  has  been 
established,  it  may  be  applied  on  a  "point-to-point  basis"  to  a  rather 
complex  ball/race  contact,  and  the  validation  of  the  model  with  the  ball 
bearing  data  shall  become  significantly  meaningful. 

10.  The  ultimate  outcome  of  any  analytical  model  of  a  mechanical  component 
is  normally  judged  in  terms  of  the  accuracy  of  overall  performance  and 
life  predictions.  With  a  validated  traction  model  the  understanding  of 
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shear  stresses,  film  thicknesses  and  the  overall  tribological 
interaction  of  mating  surfaces  in  practical  components,  such  as  rolling 
bearings,  shall  be  significantly  enhanced.  The  current  life  prediction 
models  can,  therefore,  be  reliably  updated  to  allow  for  surface- 
originated  fatigue  and  asperity  interactions  in  heavily  loaded  high¬ 
speed  bearings. 
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Appendix  A 

TRACTION  DATA 


The  data  obtained  at  the  AFWAL  Aero  Propulsion  Laboratory  and  used  in  the 
present  investigation  for  the  development  of  traction  model  coefficients  is 
documented  in  this  Appendix.  The  data  sets  are  identified  by  set  numbers  A, 
B  and  C,  for  the  different  lubricants  and  contact  configurations  as  follows: 

Set  A:  MIL-L-7808  lubricant  (Code  077-11),  Circular  Contact. 

Set  B:  MIL-L-7808  lubricant  (Code  077-11),  Elliptical  Contact. 

Set  C:  MIL-L-23699  lubricant  (Code  077-15),  Elliptical  Contact. 

The  data  is  documented  in  terms  of  selected  experimental  data  points,  indi¬ 
cated  by  various  markers,  and  model  predictions,  indicated  by  solid  lines, 
in  the  graphs  included  in  this  Appendix.  The  various  symbols  and  text  for¬ 
mat  on  each  graph  is  explained  by  the  following  labels,  indicated  in  bold 
face  on  a  typical  graph  shown  on  the  next  page. 

1.  Data  set  number  assigned  as  follows: 

First  Digit:  A,  B  or  C  indicating  lubricant  type  as  shown  above. 

Second  Digit:  1,  2,  3  or  4  Indicating  operating  temperature  number. 
Third  Digit:  1,  2  or  3  indicating  rolling  speed  number. 

2.  Lubricant  name. 

3.  Disk  radius  (inch)  in  rolling  direction. 

4.  Crown  radius  (inch). 

5.  Inlet  Temperature  (°F). 

6.  Disk  rpm  under  pure  rolling. 

7.  Model  type  1  or  2. 

8.  Reference  temperature  (°K)  used  in  the  model. 

9.  Viscosity  parameter  (Pa.S),  y0*  or  y„**. 

10.  Pressure-viscosity  coefficient  (1/Pa),  a*  or  a**. 

11.  Temperature-viscosity  coefficient,  ft*  (1/°K)  or  0**  (°K). 

12.  SSD  for  the  regression  fit. 

13.  Marker  identifying  the  experimental  data. 

14.  Solid  line  number  corresponding  to  model  predictions. 

15.  Hertz  contact  pressure  (MPa). 

16.  Contact  load  (N). 

17.  Rolling  velocity  (M/S). 

18.  Inlet  temperature  (°K). 

19.  Traction  coefficient. 

20.  Slide/Roll  ratio. 

For  data  set  A,  two  sets  of  plots  are  included.  In  the  first  set,  label  #5, 
inlet  temperature  (°F),  Includes  an  asterisk  (*),  which  means  that  the  model 
correlations  were  obtained  with  the  disk  temperature  set  equal  to  the  inlet 
temperature.  In  the  second  set,  which  is  without  the  asterisk  on  label  #5, 
the  correlations  were  obtained  with  actual  disk  temperature  data  included  in 
the  experimental  data  set.  All  correlations  for  data  sets  B  and  C  were 
obtained  with  the  disk  temperature  set  equal  to  the  inlet  temperature. 
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